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Abstract. Inclusion of spin-dependent interactions in graphene in the vicinity of the Dirac points can be 
posed in terms of non-Abelian gauge potentials. Such gauge potentials being surrogates of physical electric 
fields and material parameters, only enjoy a limited gauge freedom. A general gauge transformation thus 
in general changes the physical model. We argue that this property can be useful in connecting reference 
physical situations, such as free particle or Rashba interactions to non-trivial physical Hamiltonians with a 
new set of spin-orbit interactions, albeit constrained to being isoenergetic. We analyse different combina¬ 
tions of spin-orbit interactions in the case of monolayer graphene and show how they are related by means 
of selected non-Abelian gauge transformations. 
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1 Introduction 

The spin orbit interaction (SOI) is rapidly becoming a tool 
for spin manipulation in material systems of high techno¬ 
logical importance such as non-centrosymmetric semicon¬ 
ductors and of potential relevance such as graphene-like 
compounds [Tl[2]. It is also a crucial ingredient in the cre¬ 
ation of surface chiral transport states [31ll] of topologi¬ 
cal insulators where the bulk gap is created by the SOI. 
The interaction of relativistic origin, couples the spin to 
an electric field whose source can be internal, due to the 
atomic cores, or externally applied via voltage gates. The 
spin-orbit coupling enters the Hamiltonian in two variants, 
i) the Rashba type coupling that involves an external elec¬ 
tric field via structural inversion asymmetry[5] produced 
for example by a gate and ii) the intrinsic spin-orbit cou¬ 
pling, produced by the inversion asymmetry of the bulk 
material as in zinc-blende crystals or purely by the core 
penetration of the atomic wave functions 0. Analogous in¬ 
teractions are defined for two dimensional materials such 
as graphene, where Rashba type interactions are a result 
of a combination of the atomic SOI and the Stark effect 
due to an external field, while the intrinsic coupling is 
only a function of the lattice structure and the atomic SO 
coupling. 

The spin-orbit interaction might appear as a peculiar 
interaction to describe in a Hamiltonian formulation since 
it depends on momentum, so in a sense it is part of the 


complete description of a kinetic energy operator. In this 
sense there is a mathematical connection to gauge theory, 
where one can reinterpret the SOI as a gauge vector and 
gain insight into the problem by way of the connection to 
well established formalisms in QED and gravity [B]. 

As a direct illustration let us take the Pauli SO cou¬ 
pling where there are no external magnetic fields. 


es ■ (E X p) 

2m 2m2c2 ’ ^ ^ 

where p is the momentum operator, s is the electron spin, 
E is the electric field m is the electron mass and c is the 
velocity of light. It is immediately suggestive to complete 
the square and construct a Hamiltonian that is reminis¬ 
cent of a minimal coupling in electromagnetism. 


H = 


1 

2m 


8m 


( 2 ) 


where = {e/2mc^)eajiEj can be seen as a new SU{2) 
gauge field and the definition of a minimal coupling to the 
spin is done by the substitution Pi ^ Pi — This 

is equivalent to defining the covariant derivative = 
di — Here, the SU{2) gauge field is IT“ and the 

generators of 51/(2) algebra, given in terms of the spin 
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components Sa where Sa = ha a/‘2- and a a is a Pauli ma¬ 
trix. In the expression W^Sa^ the spin index a is contracted 
and the matrix structure of the non-Abelian gauge field is 
apparent. The space index i runs over i = 1, 2 in standard 
applications to two-dimensional non-relativistic electron 
gas (2DEG). It is important to note that in this formula¬ 
tion the gauge field is anchored to physical quantities such 
as material parameters and the electric field [5]. 

Such a gauge point of view, reformulating in general 
terms the spin-orbit interaction as a SU{2) x ?7(1) gauge 
theory in the non-relativistic case, has been known for 
some time [711H1IS]. Specific applications in the case of the 
Pauli equation and the Rashba and/or Dresselhaus spin- 
orbit interactions (SOI) in 2DEGs were analysed in refs. [101 
[TTlll2U13l[T^[T51[TMr7] . Applications of SU{2) gauge for¬ 
malism to the case of transport in a superconductor with 
SO coupling is studied via Green functions formalism in |18j . 

Because the non-Abelian gauge fields components IE“ 
are proportional to the electric field, gauge transforma¬ 
tions U of the system’s Hamiltonian H amount to trans¬ 
formations to a form UHU~^ that is unitarily equivalent, 
hence iso-spectral, but describes in general a different dy¬ 
namical situation. Indeed, non-Abelian gauge transforma¬ 
tions IE“ —>■ IE“ -I- — Cabco’^W^ contain a term corre¬ 

sponding to an internal rotation. In the case of the SO 
gauge field, the internal space is spin space and a ro¬ 
tation in spin space as described by the term tabcO^Wf 
in the gauge field transformation is precisely what SOI 
does: it produces spin precession and is explicitly used 
to manipulate spin orientation in 2DEGs. In the case of 
Abelian theories, one needs singular gauge transforma¬ 
tions to produce a new physical situation by a modifi¬ 
cation of the curvature m, here it may be realized even 
with a global gauge transformation. Nevertheless, there 
remain classes of gauge transformations that preserve the 
physics, that have been exploited in several papers in 
the context of 2DEGs [20lini[2^[2m[24l[^[26] . In the first 
of these papers m a two-arms spin interferometer with 
Rashba SOI is reduced to a system with two indepen¬ 
dent C/(I) fields through a SU{2) gauge transformation. 
In ref. [H], a 2DEG with both Rashba and Dresselhaus 
SO interactions is studied in the Abelian limit of equal 
SOI amplitudes. Chen and Chang then use an Abelian 
gauge transformation to map the system on a free parti¬ 
cle model and they recover the phenomenon of Persistent 
Spin Helix (PSH) first analysed by Bernevig, Orenstein 
and Zhang [3]. Other authors have already used gauge 
transformations to simplify SO-related problems. Levitov 
and Rashba consider SO coupling in an anisotropic semi¬ 
conductor quantum dot where the behaviour of electrons 
is governed by a Schrdinger equation. A gauge transfor¬ 
mation was used to turn the SOI into a Zeeman interac¬ 
tion [23]. In [23], Tokatly and Sherman consider a 2DEG 
with SO interaction where they analyse diffusive spin dy¬ 
namics and precessional spin dynamics. A gauge transfor¬ 
mation is used to gauge away a pure gauge SU{2) field. 
In [21] the authors study a non relativistic 2DEG with 
weak inhomogeneous SO coupling to propose a spin tran¬ 


sistor which works not only in the ballistic regime, but 
also in the diffusive regime. Their approach is based on 
Onsager reciprocity relations, and they use a local SU{2) 
gauge transformation to get a spin diagonal structure in 
the perturbative regime. In a similar perspective in |26j 
the authors analyse a 2DEG with SO coupling in relation 
with Onsager relations to relate the vanishing of spin Hall 
conductivity when there is only Rashba SO and the finite 
inverse spin Hall effect. An SU{2) gauge transformation 
allows them to shift the SO coupling to a Zeeman inter¬ 
action. 

A similar gauge formulation is valid in the relativis¬ 
tic case leading to a theory with a structure very simi¬ 
lar to the electroweak theory m ( in two space dimen¬ 
sions [2H]) where C = Here, S' is a two com¬ 

ponent spinor. The Dirac matrices satisfying Clifford 
algebra { 7 '", 7 ‘'} = 2 ? 7'^'^1 (with = Diag[l,-1,-1]), 
are at least 2 x 2 matrices [22] and a representation can 
be given in terms of Pauli matrices 7 ® = Uz, 7 ^ = iay, 
7 ^ = —ia^. The Dirac equation in this case is a 2 x 2 
matrix equation which, in the absence of an external U{1) 
field, takes the form 

Hd = ca. ■ p + l3mc^, (3) 

where one identifies = a^ and = ay, and j3 = a^ 
( 7 ° = /3 and 7 = (3a) . Like in the non-relativistic case, the 
non-Abelian gauge field components will follow from the 
identification of the interactions coupled to the electrons 
spins through the minimal coupling. 

In contrast to the non-relativistic case, the dispersion 
relation of massless Dirac particles is linear in momen¬ 
tum so the SOI, can be embedded as a new SU{2) gauge 
field. As the new concocted gauge fields, built by anal¬ 
ogy to a gauge formulation are also fixed by the phys¬ 
ical situation, we can then exploit the whole machinery 
of gauge transformations in a completely different fash¬ 
ion (not contemplated even for quadratic dispersions): In¬ 
stead of seeking special classes of transformations that pre¬ 
serve the physics, we determine the gauge transformations 
that change it (because material parameters and fields get 
modified), in order to consider new non-trivial physical sit¬ 
uations which, to our knowledge, have not been addressed 
directly. This is the case of an in-plane electric field (the 
Rashba interaction corresponds to an out-of-plane electric 
field), or the case of non uniform (i.e. spatially dependent) 
gauge parameters. 

The paper is organised as follows. In section 2 we re¬ 
visit the Hamiltonian of single layer graphene with SOI 
and show how it can be rewritten in terms of a fictitious 
gauge field. In section 3 we choose a class of gauge trans¬ 
formations that mute the regular Rashba SOI auxiliary 
gauge field into a tunable combination of Rashba, Dres¬ 
selhaus and in plane electric field contributions. Section 
4 is devoted to treat the problem of connecting the free 
graphene without interactions into a spatially dependent 
SOL Finally, section 5 addresses the generation of Stark 
interaction in graphene [30] . As these transformations are 
done exactly we also obtain the wave functions for a com¬ 
pletely new physical situation, thus completely solving the 
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problem with new couplings. The last section summarizes 
the paper. 

2 The graphene Hamilonian 

The Kane-Mele Hamiltonian [3] , a continuum Hamiltonian 
derived from a tight binding formulation expanded around 
the K points, includes both intrinsic and Rashba spin- 
orbit interactions, 

Hk = uf ct • p -b Asoo-zSz + >^K{cr^Sy - ayS^). (4) 

uf is the Fermi velocity, the parameter /Iso measures the 
strength of the intrinsic SO coupling, and Xr is the Rashba 
spin-orbit strength whose source is a Structural Inversion 
Asymmetry (SIA) and is proportional to an external out- 
of-plane electric field. The origin of term springs from a 
tight-binding approach that couples nearest-neighbour pz 
orbitals through paths that involve Pz — Px,y coupling to 
the graphene plane, a Vgpa hopping integral, and then a 
Stark s—pz coupling. So the Xr involves both the external 
electric field and the material parameters. 

In order to introduce gauge fields, we use the mini¬ 
mal coupling and write equation (| 3 ]) in the form of uf o’ • 
(p — W“ So)- There is nevertheless a caveat there, since 
this expression is not Hermitian if a does not commute 
with p — W“ Sq. The recipe in that case is to work with 
symmetrized expressions. Once properly symmetrized to 
generally account for coordinate transformations |15) , the 
Hamiltonian can be written as 

Hk = VF (T ■ {v- W“ Sa) -b iuF[p - W“sQ,tT], (5) 


an interaction already present. We can transform the dy¬ 
namical momentum tt = p — W“sa which follows from 
the minimal coupling prescription, under a non-Abelian 
gauge transformation parameterized by the set of three 
functions (see e.g. Ref. [32]) in Us = 6^^®“ according 
to 

p - W'%, =w' = Us TT ur^ 

= p _ (W“ + - eabca^W^) Sa (7) 

up to 0(|ap) terms, where |a| = y/a’^Ua- The trans¬ 
formation of the gauge potential thus assumes the form 
W'“ = W“ -b VQf“ — eabcO^W'^. Even a set of constant 
gauge transformation parameters a“ will in general create 
new interactions W'“ from original ones W“ through the 
last term of the equation. On the other hand, if one starts 
at the outset from the free particle case, non uniform Q;“(r) 
are required to generate non trivial effects. 

Let us consider a graphene sheet and contemplate an 
initial problem with a non zero spin-orbit interaction, e.g. 
Rashba SOI due to the presence of an electric held E± 
perpendicular to the sheet (here treated in Cartesian co¬ 
ordinates at one of the Dirac points), 

Hr = vf o--p + XB.{axSy-aySx). ( 8 ) 

Acting on a spinor in the form of a two-dimensional plane 
wave multiplied by 4-component space-independent am¬ 
plitudes, the Hamiltonian of Eq. (I8|) leads to the secular 
equation {E/hvp)'^ — ^{E/hvF)'^l\^\'^ +2 {Xr/vf)'^) + \^\'^ = 
0. The eigenvalues follow m 


where the last commutator is —ih (V - a). When both Rashba 
and intrinsic SOI are present, inspection of Eq. o shows 
that the gauge helds take the form 

= (Ar/uf)u^, W2 = -(Ar/uf)u„ 

W3 = -{Aso/vf)Uz, (6) 

where are unit vector in the f-th cartesian coordinate. 

Ar and Aso both have the dimensions of inverse time and 
IF“ those of an inverse length. 

We now reemphasise that the present gauge fields de¬ 
pend on physical given quantities so that there is no gauge 
freedom [13U31j (an introductory review can be found in 
Ref. [32]) ■ We take advantage of this peculiarity of the 
gauge structure of the spin-orbit interaction to find a gauge 
transformation that connects two physically distinct prob¬ 
lems. This we develop in the next section. 


Ek.,s = nhvF (^SXr/vf + \/|kp -b (Ar/uf)^^ , (9) 

where k = ±1 is the particle-hole index and ^ = ±1 labels 
the SO branches. 

We will investigate the transformation Ug = exp 
of this Hamiltonian under a set of constant gauge param¬ 
eters a“. In this case the gauge operations are pure ro¬ 
tations of the spin degree of freedom so one can generate 
very simple interaction transformations. The free particle 
term is invariant and we have to consider the transforma¬ 
tion of terms like 

i q;“ A ( i \ 
--s.jsbexp(^---s,j 

2 |ct| ^ |a| • |a| a 

= COS Vsb - — cos Vsm yeobcQ! Sc 

|a| ^ ^ 


3 Transforming spin-orbit couplings on a 
graphene sheet 


(a)^ 


■ sin^ a°^a'^{5bcSa - SacSb + SabSc)- 


( 10 ) 


We will now contemplate the case of spin-orbit interaction 
which can be described in terms of non-Abelian gauge 
theory. Our aim is to generate an interaction or modify 


A case which has enough generality for our purpose is 
given via the choice = a cost/), o? = asin^, o? = 0 , 
(/> = const. A rich variety of spin-orbit interactions are 
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generated, e.g. when we fix </> = 7r/4, 


one obtains the eigenstates of Hamiltonian (na). 


Ar(' 


(JxSy CT-J/S 


y^x 


) -f Ar COS^ ^(CT; 




x^y 


y^x 


-I-Ar sin' 




A 


Sy) H —ji sin a(cTa 
v2 


( 11 ) 


^\V) = 


I 


= e 


ikr 


ie cos ■ 


cos ■ 


H 

2 

k (5 cos 1 ^ sin 


■z^sin 


sin ■ 






_-i^g-i0k sinMcos%- 


^ V 


|a:| 

I COS ^ sin 


+ cos 


2 

^ n d I 

cos^J 


i.e. one generates a Dresselhaus interaction with ampli¬ 
tude Ap = ARsin^a/2 and a Pauli spin-orbit interac¬ 
tion which comprises a term due to a perpendicular elec¬ 
tric field (ordinary Rashba term with amplitude A(^ = 

Ar cos^ ci/2) plus another term due to an in plane electric 
held with amplitude Ay = ^ sin a. 

Interesting limiting cases are found when the ampli¬ 
tude a is hxed: 

i) a = TT, which amounts to a rotation by this angle 
around the x = y axis, where the transformed problem 
exhibits a pure Dresselhaus interaction, and 

ii) a = 7r/2 where Rashba and Dresselhaus have equal 
amplitudes and the in-plane Pauli interaction is also present 
in the transformed problem. Case i) interestingly shows 
that Rashba and Dresselhaus interactions are related through 
a gauge transformation. This case has been addressed in 
the literature [34]. 

For case ii) the transformed Hamiltonian reads as 


(15) 


F^R — Up O’ • p -|- Aj^(tTxSy — CTySx) 




O', 


iSy') XtAcfxSz ^y^z) 


The eigenvalues are the same as those of EqlHl As the phys¬ 
ical fields and the material parameters generating the cou¬ 
plings in the above Hamiltonian are different from those 
of the original problem, we have then fully solved a new 
physical situation. 

There is nevertheless a physical effect that we have 
not taken into account. When there is an in plane electric 
field (as the one which enters the spin-orbit interaction 
in the transformed model), one has to add the Stark con¬ 
tribution to the potential energy, which, although spin- 
independent, has a spatial dependence which modifies the 
physical problem under study. This example allows us to 
illustrate another route to generate spin-orbit interactions 
from the gauge transformation. 


4 Spatially dependent transformations 


We now contemplate a graphene sheet and start from the 
free particle problem (at one Dirac point) in Cartesian 
coordinates. 


( 12 ) 


Hq = vp a ■ p = -ihvp{axdx + cFydy), 


(16) 


The corresponding (un-normalized) eigenvectors of the ini¬ 
tial Hamiltonian are given by 


iF(r) = 


CQg 

k5 sin 
i sin 

cos / 


(13) 


with tan())k = ky/kx and tanS^.i = EK,,s/hvp\\i\, 0 k,, s € 

[-7r,7r]. 

Acting on these eigenstates, with the unitary gauge 
transformation 


D, = 


/ 

|q:| 

COS ^ 

0 

^sinM 

0 


0 

|a| 

COS ^ 

0 




0 

|a| 

COS ^ 

0 

V 

0 


0 

|a| 

COS ^ 


(14) 


and assume a gauge transformation Ug = exp(j-^Sa) 
with a space-dependent set of Q;“(r). The idea is here to 
generate local SOI by rendering local the gauge structure 
of the theory. The gauge transformation of the first term 
in Eq. m reads as 


(idx) ^ exp{^^Sa)iidx)exp{-^Y^c)- (17) 

The transformed operator at the r.h.s. of Eq. being a 
2x2 matrix, it can be rewritten in terms of Pauli matrices 
(in real spin space): 


idx idx + 

1 . 


1 


a 

-sm —- 

a 2 


(a^a") - 


a^(dxa) 


a 


W Sc . |a| _ ^ Sb. .. 

COS-^— - sm^(^acls -^eabc^) ^^) 

2 a a 2 a 


This expression is very involved, but it clearly shows how 
diverse forms of spin-orbit interactions may be realized af¬ 
ter the tensor product between the pseudo-spin matrices. 
As an example, let us consider a space-dependent gauge 
parameter consisting in only one spin space component 
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which depends on the position in the graphene sheet, 
= a^{x, y). Then, the free particle Hamiltonian vp (^ P 
in Eq. m becomes 

Hq = vp a ■ p - + dya^UyS^) 

= vp a ■ {p-'W°‘{x,y)sa). (19) 

The complete solution for the “in-plane field” Hamilto¬ 
nian da can now be obtained via the unitary spatially 

dependent gauge transformation [/s = exp(|- ‘^ ^2 ’^^ Sz) ap¬ 
plied to the free problem (ITSl) . An example is treated in 
the next section. 


5 The Stark coupling 

The transformation performed in section [3l generates an 
in-plane electric field £ = {£x,£y)- Then, the associated 
contribution to the potential energy must be taken into 
account. For simplicity, we consider the case of a uniform 
electric field. The Stark potential energy —e£ ■ r, in terms 
of a tight-binding Hamiltonian [35], and then expanded in 
the vicinity of the Dirac point in the continuum limit, 
reads as Hstark — ^xx^xPx ^xy^xPy ^yx^yPx ^yy^yPy 
(up to a quantity proportional to the identity that we can 
absorb in the redefinition of the zero of the energies). The 
coefficients are given in terms of the hopping elements 
along the bonds of the lattice that are now anisotropic and 
are labelled as 7^: ujxx = ^'^(72 +73), ^^xy = ^'^(72 - 

73), U}yx = -^(72-73), Wyy = ” f ( 7 i T ^ (72 T 73 ))■ Ooly 
the form of the Hamiltonian is important here. 

When the Stark contribution is added to the ordi¬ 
nary kinetic energy, one is thus led to an anisotropic ki¬ 
netic energy which reflects the drift induced by the field: 
‘^xx^xPx '^xy^xPy ‘^yx^yPx '^yy^yPy with = Vp -\-UJii 
and Vij = ujij. An interesting simplification occurs when 
the in-plane electric field is along the 71 bond, since then 
72 = 73 and Wxy = ujyx = 0. If now one considers the 
consistent Hamiltonian comprising the ordinary kinetic 
energy and the Stark and spin-orbit contributions asso¬ 
ciated to an in-plane electric field £ = £ uy , one has 

CVp 

H = Vxx ^xPx T '^yy ^yPy “f ^ (^9) 


The cancellation of the spin-dependent terms is achieved 
by choosing a^(r) = — 2 m Jv ^ and the Hamiltonian thus 
becomes 

H - Vxx ^xPx Vyy (JyPy. (^^) 

Its eigenvalues and eigenvectors are easily obtained. 


_ ^i(k^x+kyy) 


/F^ 

n 




(23) 


E' = ±Y^ {hvxxkxY + {f^VyykyY, (24) 

and those of the original problem (EOj) then follow via the 
inverse gauge transformation 


xpiv) = 


e£vp 


Amc'^vx 

e£vp 

“ 

Arnc'^vxx 


F^e 

G-^e 

\Gie 


A ^ 

<^mc^Vxx 
_. eSvp 


■ (25) 


The physical properties being associated to different 
physical situations, they may differ in the two gauges. 
Consider the example of the current densities. The charge 
current density in a given eigenstate \F) may be defined 
by the matrix element jq = e{F\vp a\F) where vp a = 
{vp)~^{vxxO'xVix + VyyayUy) Is the velocity operator. This 
charge current density is invariant under the gauge trans¬ 
formation, since in the primed problem one has 

j'q = evp{F\U~^{axUx + ayUy)Us\F) =jq (26) 

because the operator in sandwich is spin-independent. 
This is not a surprising result, since as we have noticed 
the non-Abelian gauge transformation produces a spin ro¬ 
tation and hence does not touch physical properties which 
are independent of the spin state. The situation is differ¬ 
ent with the spin current density, j“ = {F\vp aSa\F), since 
then Ug does not commute with Sa and the resulting spin 
current is depending on the gauge. Like in the case of the 
charge current, let us form the spin current in the primed 
problem in terms of the original eigenstates, 

= {F\U-^vp asaUg\F) (27) 


Devising a transformation that eliminates the space-de- 
pendent SOI term, we look for one that is the inverse to the 
one that created the interaction i.e. Ug = exp . 

Using Eq. m, one gets Px,y -)► Px,y - ^{dx,ya^)sz and 
the transformed Hamiltonian is now 


H - UgHU^ - Vxx ^xPx Vyy (7yPy 

'1 


3 evp£ \ 
XXGfx^ “T . o I 
4??1C^ J 


2 
1 

-VyyWy 


^X^Z 


- I ly d 


(X ) (JySz' 


( 21 ) 


where Ug ^SaUg -)■ cos^ ^Sa + sin^ ^(2s^5a3 - Sa), i-e. 
Sz Sz because the gauge transformation consists in a 
rotation around the z—axis, while Sx,y —>■ cos^ — 

sin^ An example of a spin-dependent physical prop¬ 

erty which is not the same in the two problems considered 
is 

'fg = {F\vpaSx\F) (28) 

with Hamiltonian (l20l) and 


{F\vp a 



\F) (29) 
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with Hamiltonian (1221) . 

This is a known result in non-Abelian gauge theory 
that the conserved currents are gauge dependent (gauge 
covariant) [dfi] . The spin current furthermore constitues 
only part of the conserved current and, alone, is not con¬ 
served due to the fact that the spin interacts with the 
non-Abelian gauge fields (i.e. here the electric field), see 
e.g. Ref. [37]. 


6 Summary and Conclusions 

We have analysed SU (2) gauge transformations in graphene 
close to the K points. As unitary transformations they 
preserve the eigenvalues, but due to the electric field de¬ 
pendence in the auxiliary gauge potentials, such trans¬ 
formations convert from one physical model to another. 
Here we contemplated the connection between the pure 
Rashba coupling to the pure Dresselhaus coupling and 
also between the former and a Dresselhaus coupling plus 
an in-plane electric held. We also showed how one could 
generate SOI resulting from in-plane electric helds from 
the free graphene Hamiltonian thru a spatially dependent 
gauge transformation. The strengths of the couplings were 
proportional to the gradient of the space dependent gauge 
parameter. Physically, the in-plane Held effect can also be 
realised in practice for stretched graphene. 

We Hnally considered a series of tailored gauge trans¬ 
formations in order to address the consistent Hamiltonian 
for an in-plane Held plus its corresponding Stark cou¬ 
pling. in this case, we transformed the problem into an 
anisotropic Dirac cone problem with intrinsic SOI. For 
this problem, the eigenfunctions can be fully derived and 
transformed back to obtain the corresponding ones for the 
original Hamiltonian. 

We must note that the gauge transformation fixes the 
relative values of the interactions generated as a function 
of the gauge transformation chosen and the original pa¬ 
rameters. If there is a threshold for a phase transition as a 
function of the parameter values, the transformation will 
in general not address it. 

We have stressed that the SU{2) gauge transforma¬ 
tion applied to the Dirac Hamiltonian, for the three prob¬ 
lems considered before, now corresponds to different phys¬ 
ical situations, since the non-Abelian gauge field compo¬ 
nents are bound to real external fields and material pa¬ 
rameters. In this manner, we have shown that the cor¬ 
responding Hamiltonian can be conveniently addressed 
by means of the aforementioned applied gauge transfor¬ 
mation. The transformation however keeps the eigenener- 
gies unchanged, so the modification of eigenvectors and of 
spin-dependent operators after gauge transformation now 
allows for different physical expectation values regarding 
for example spin persistent currents. This is connected to 
the covariance property of SU{2) gauge transformations 
which amount to rotations in internal space which here is 
the ordinary space coupled to spin degrees of freedom. 
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